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We investigate low-energy dynamical properties of one-dimensional multicomponent quantum 
liquids with the short-range interaction as well as the l/a;-type long-range interaction. By 
calculating the single-particle spectrum and the dynamical spin susceptibility by means of 
the bosonization method, we discuss how the orbital degeneracy and the band splitting affect 
the dynamical response functions. The effect of the long-range interaction is also addressed. 
Although the long-range interaction suppresses charge fluctuations, it effectively enhances spin 
fluctuations via the formation of the Wigner crystal. 
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§1. Introduction 

Low-energy excitations in one-dimensional (ID) correlated electron systems with the short-range 
interaction are described by collective charge and spin fluctuations. This class of quantum liquids is 
referred to as the Tomonaga-Luttinger (TL) liquid.Bi'i-* Recently, dynamical properties of ID elec- 
tron systems have been studied extensively by various experimental methods, e.g. angle resolved 
photoemission spectroscopy (ARPES),B& NMR|li'& etc. In particular, the spin-charge separa- 
tion inherent in ID electron systems was observed in the single-particle spectrum via the ARPES 
measurementsj^i^ and typical power-law behaviors in spin correlation functions were observed in 
the temperature dependence of the NMR relaxation ratei'BIHI' in various ID compounds.i'0i) 

More recently, the effect of the orbital degrees of freedom has attracted renewed interest in cor- 
related electron systems. Such orbital effects are indeed important in ID. For example, in ID 
correlated electron systems at quarter-filling, such as Na2Ti2Sb2O0 and NaV2O5,lll0 
the effective orbital degrees of freedom are essential to understand the correct low-temperature 
properties. In this connection, the ID multicomponent spin-orbital model has been extensively 



* E-mail : satoshi@tp.ap.eng.osaka-u.ac.jp 



studied analyticallySSiliiili'i) and numerically. 011) However, dynamical correlation func- 
tions have not been studied systematically so far, except for exactly solvable 1/r^-models. 

It is usually assumed that the Coulomb interaction is well screened, and the resulting interaction 
is short-ranged in the TL liquid. However, this is not always the case for quasi-lD systems, since 
the screening effect due to adjacent chains is not efficient enough in some cases, leading to the bare 
long-range interaction. The resulting quantum liquid is a ID analog of the Wigner crystal, which 
exhibits quite different properties from those of the TL liquid. 

In this paper, we investigate the effect of the orbital degeneracy on low-energy dynamical prop- 
erties of the TL liquid by employing the bosonization method for the two-orbital electron model. 
We also clarify how the long-range interaction affects dynamical properties via the formation of the 
ID Wigner crystal. This paper is organized as follows. In §2, we introduce the model and derive a 
low-energy effective Hamiltonian in the bosonized form. In §3, we investigate the orbital effects on 
the TL liquid with the short-range interaction and then in §4 discuss the effect of the long-range 
interaction on the low-energy properties of the system. A brief summary is given in the last section. 

§2. Low-Energy Effective Hamiltonian 

In this section, we outline the derivation of a low-energy effective Hamiltonian in the bosonized 
form. We start with the ID two-orbital Hubbard model with the band splitting A, 



n = - t^{c]^^Cj+laa + h.c.) 



i,j a,b a, a' 

- [(njii + rijii) - K-at + nj2i)], (1) 

where c^^^^ creates an electron with the orbital index a = 1 (lower band), 2(upper band) and spin 
o" = 1, 1 at the j-th site. In the absence of the band splitting, A=0, the above model corresponds 
to the SU(4) symmetric Hubbard model. For simplicity, we have neglected the Hund-rule coupling 
between different orbitals, which are irrelevant in the following discussions in the TL liquid phase. 

Passing to the continuum limit, we derive a low-energy effective Hamiltonian in terms of the 
slowly varying right- and left-going Fermi fields Cj^l- and CjJ^, which are related to the original 
lattice operator as 

Cjacr ~ V"o e ^ Cj^^ -h e Cj^^ , [Z) 

where oq is the lattice spacing and kp^ (kp^) is the Fermi momentum of the lower (upper) band. 
Since we are interested in a metallic phase with the repulsive interaction, we are left with only 
forward scatterings. By using the Fourier decompositions Cj"]^ = -^J2p^~^^''"'°Crpaa {r=+ or -), 
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the Hamiltonian (1) is rewritten as 



p,a,a 
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a.a' 



(3) 



where Praa{p) = J2k ^Ik+paa^rkaa are the density operators satisfying the bosonic commutation 
relation, and Vp^ = agt sm{kp^ gq) is the Fermi velocity of the lower band (a = 1) or the upper 
band (a = 2). We have four kinds of bosonic fields, pri-[, Prij.! Pr2'\ smd Pr2j,i which are converted 
via a unitary transformation to the new bosonic fields corresponding to the 'charge' (c), 'spin'(s), 
'orbital' (o) and 'spin-orbital' (so) sectors, 
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(4) 



Moreover, when the band splitting takes a finite value, it is convenient to use the new basis for the 
charge and orbital sectors, 

p+c + P~c I ^ / ^^(P) -^3{p) \ I p+c + P-c 
p+o + P~o J \ Clip) J \ p+o + P-o 

where 



(5) 



= cos[a(p)]. 



(6) 
(7) 
(8) 



The explicit formulae for y{p) and a(p) are given in Appendix. Note that when two orbitals are 
degenerate (A=0), we have a{p) = 0, so that = 1 and ^2 = ■^3 = 0. 
Consequently, we end up with the Hamiltonian, 



^ = E E ^-(P) I -t;^[p+Ap) + P-Ap)? +K,{p)[p+,ip) - p-uip)?} , 



(9) 



where the velocities and the TL parameters arc summarized in Appendix. 

This completes the derivation of a low-energy effective theory of ID correlated electrons with two 
orbitals. Though we have outlined the way in terms of the Hubbard model, the resulting formula 
can be applied to generic two-band models in a metallic phase. 
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Before closing this section, some comments are in order for the formula (9). In terms of the spin 
and spin-orbital basis, the velocities read Vs=Vp^ and Vso=Vp^ ■ Moreover, we should set Ks=Kso=^ 
according to SU(2) symmetry required for these excitation modes. 

§3. Multicomponent Tomonaga-Luttinger Liquid 

We first consider the system with the short-range interaction and discuss how the orbital degrees 
of freedom affect the dynamical response functions. When the interaction is short-ranged {6- 
function), we set all the parameters, Viy{p) and K^{p), independent of the momentum p. We 
employ the i- J model as a microscopic model to give the values of v'^ ^ and K[. ^. Furthermore, if we 
make use of the integrable supersymmetric i-J model,lii'0'ii'0^ where the hopping t and the spin 
exchange coupling J satisfies the relation, t = J, the above TL parameters can be determined by the 
Bethe ansatz method combined with conformal field theory. This model is somewhat special, but 
it can well describe general properties of the doped spin-orbital model, as far as we are concerned 
with a metallic system near the insulating phase. 

We first derive the single-particle Green functionBSiiS) for the TL liquid with orbital degen- 
eracy, which is defined as 

GraAx,t) = -ie(t)([Vraa(x,t),VL(0,0)] + ) 

= -ie(t)[G<,(x,t) + G<,(-X,-t)], (10) 

where 

GraAx,t) =< V'raa(x,t)4aa(0,0) >, (11) 

and G^^^(x,t) is defined similarly. The Fermi field is written down as 

^^'^'^ = i^o (-^ ? - ^''] ) ' ^^^^ 

where ijraa is the Klein factor, which ensures the anti-commutation relation between the right- and 
left-going electrons. The above two phase fields are 



L '-^ p 



where 



eL.(p) = 6(p) + (-i)"+'e3(p) 



(13) 
(14) 

(15) 
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a.(p) = -6(p) + (-ir+^ei(p) ' 
cra.b) = eiL(p) = ^{i-(-ir+H ■ 

Note that Oy is the bosonic field conjugate to (f)^, which satisfies [(j)y{x), 9yi{x'y\ = \'^5yysgn{x — x') . 
Inserting these expressions into the Green function, we have 
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G<.(x,t) = - n 



u=c,s. 
o,so 



sinh 



-if 2A+sgn(x+i}^t) if 2A-sgn(x-5„t) / tt V 



sinh 



2A: 



(18) 



with P = 1/kBT, where the conformal dimensions are given by 
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2A± = -[K^C^J + Cea.VK,T2re4,aM,J, (19) 

and the corresponding scahng dimensions are Xjy=A++A^. For A=0, we have C'^aa~^ea(7 which 
gives the conformal dimensions for the SU(4) spin model0 Equation (18) implies that the spin- 
charge-orbital separation occurs in elementary excitations, as should be expected for ID electron 
systems. Note, however, that all these degrees of freedom are combined together to contribute to 
the single-particle Green function. Fourier transformation of (10) gives the single-particle spectrum, 

G'r-ao-(^) <^)- 

Let us now investigate the single-particle spectrum by focusing on right-moving fermions (r=+). 
To see the effect of the orbital degeneracy, we first compare the spectral function of the single- 
orbital electron system (2-component TL liquid) with that of the two-orbital one (4-component TL 
liquid). The system has the dispersions, uj=Vuk, so that the spectral function may have several 
peaks around uj=Vyk at finite temperatures. We use the parameters which are derived from the 
exact solution of the multicomponent supersymmetric t-J modeEi) near the Mott insulator in 
order to parameterize the velocities and the TL parameters as a function of the band splitting. For 
the case of single- and two-orbital TL liquids, we set Vs=2TT/n, K^. = 1/n (n = 2 or 4) whereas 
Ks=l is fixed due to SU(2) or SU(4) symmetry. We set a small value for the charge velocity 
(t;c=0.1) because we are interested in the vicinity of the Mott insulator. Note that we will measure 
the energy in the unit of the spin exchange coupling J. 

We show the results for the TL liquid with and without orbital degeneracy in Fig. 1. At high 
temperatures, all the kinds of excitations show over-damped behaviors in which the spectral weight 
is distributed broadly around w ~ 0. This is valid both for the single- and two-orbital models. 
As the temperature decreases, the peak structure of the spectrum is gradually developed around 
each mode of elementary excitations. In the single-orbital case, the peak around the charge mode 
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Fig. 1. Single-particle spectrum for the TL liquid; (a) SU(2) single-orbital case, Vc—Q.l, Vs=n, (b) SU(4) two-orbital 
case, Vc—0.1, Vs=Vo= 5so=7r/2. 



(a;~0.01) is enhanced more prominently than that for the spin mode (l(J~0.3) at low temperatures. 
On the other hand, in the two-orbital case, the peak structure for the charge mode (a;~0.01) is 
rather suppressed, while that for the spin mode (u;~0.17) is enhanced. This is understood by 
looking at the zero-temperature behavior of the spectrum. At T = 0, the above spectral function 
shows power-law divergence near each excitation mode as, 

ImG(A;,u;) ~ Iw -v^^r^^^, w > 0, (20) 

for the mode specified by where the corresponding critical exponents are listed in Table |^. It is 
seen that the spin exponent 2A~ is increased in the two-orbital case due to the orbital degeneracy, 
resulting in the enhanced peak structure. On the other hand, the charge exponent 2A~ is decreased 
in the orbitally degenerate case, being consistent with the suppressed charge mode in Fig. 1(b). In 
this way, the orbital degeneracy affects not only the internal spin- and orbital-degrees of freedom 
but also the charge degrees of freedom, changing the structure of the spectral function rather 
considerably. 
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We now discuss the effect of the band sphtting. The single-particle spectral function is shown for 
the band splitting A~0.6 in Fig. 2, which should be compared with the case of A=0 in Fig. 1(b). 
When the band splitting is present, the velocities take different values. We set these parameters as, 
Us~2.6, Wo~0.8, Uso— 0.25, Kc = 1/n, i^o=2/(nx 0.436) (n = 4) according to the exact solution of 
the supersymmetric J model near the Mott insulatorii'0) and Ks=Kso=^ due to the symmetry 
requirement. In order to clarify the properties of the spectrum, we glance at the critical exponents 



Table I. Critical exponent of the single-particle spectrum. 
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again at zero temperature listed in Table ^. It is seen in Fig. 2 that the introduction of the band 
splitting alters the shape of the spectrum rather dramatically. Firstly, we notice that the peak 
structure of the charge sector (ci;~0.01) is suppressed (enhanced) for the lower (upper) band in the 
presence of the band splitting. This is indeed seen from the change in the critical exponents; 2A~ 
decreases (increases) for the lower (upper) band with the increase of the band splitting, which has a 
tendency to smear (pronounce) the singularity around uo = Vck at T = 0. The remaining two peaks 
in the spectrum of the lower band in Fig. 2(b) are respectively identified with the contribution 
from the orbital (lj~0.08) and spin (a;~0.26) degrees of freedom. Note that the spin-orbital degrees 
of freedom do not contribute to the spectrum of the lower band. On the other hand, the spectrum 
shows a quite different behavior for the upper band (Fig. 2(a)). Namely, it has two peaks, which 
are mainly contributed by the charge and the spin-orbital (w~0.025) sectors. As seen from the 
Table ^, the orbital sector also has a finite contribution, but its critical exponent is too small to 
make a visible peak structure around T = 1/10. In this way, the introduction of the band splitting 
brings about a considerable change in the profile of the spectral function, which is quite contrasted 
to that expected for the Fermi liquid system. 
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Fig. 2. Single-particle spectrum for the two-orbital TL liquid with the band splitting A~0.6; Vc=0.1, Vs—2.6, Vo—0.8, 
Vso— 0.25, (a)upper band, (b)lower band. 



We now turn to the dynamical spin susceptibility. We show the results calculated for lm[x{k,uj)]/LO 
at and near the Mott insulator in Fig. 3. We also compare the results for the single- and two- 
orbital TL liquid (A=0). As is clearly seen from the single-orbital case, a considerable amount of 
the spectral weight is transferred from the spin part to the charge part upon hole doping. This 
dramatic change is in contrast to our naive expectation that the charge excitation has little effect 
on the dynamical spin susceptibility. It is seen from Fig. 3(b) that this tendency is also the case 
for the two-orbital model, although the rate of the transferred weight is somehow small. These 
remarkable results imply that the NMR relaxation ratei'S'Sfli'l) may be possibly enhanced upon 
hole doping via a relaxation process mediated indirectly by low-energy charge channels. In Fig. 4 
the effect of the band splitting is shown (A~0.6). In this case, the spectral weight is transferred to 
low-energy charge and orbital modes upon hole doping. 

We have so far seen that hole-doping induces the considerable transfer of the spectral weight to 
the low-energy charge mode in the dynamical spin susceptibility. In the following section, however. 
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this statement is shown to be valid only for the case of the short-range interaction, and a completely 
different behavior shows up for the long-range interaction model. 




Fig. 3. Dynamical spin susceptibility at (left) and near (right) the Mott insulator; (a) SU(2) single-orbital case, 
Vc=0.1, Vs=T^, (b) SU(4) two-orbital case, iic=0.1, Vs=Vo= Vso='^/'2. The momentum k measures the difference 
from 7r/2. 



§4. Multicomponent Wigner Crystal 

Let us now discuss the model with the long-range interaction, which may be important for 
quantum wires as well as carbon nanotubcs. Also, for ID highly correlated metallic systems close 
to the Mott insulator, the screening effect of the bare interaction may be rather poor, so that the 
long-range interaction is expected to play an important role. We will mainly focus on the latter 
case in the following discussions. 

It is known that physical properties in the presence of the long-range interaction are quite different 
from those of the TL liquid. For example, the charge density correlation function has a power-law 
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Fig. 4. Dynamical spin susceptibility for the two-orbital TL liquid in a metallic phase with the band splitting 
A ~0.6; Sc=0.1, Vs~2.6, 5o=~0.8, Vso^O.25. 

decay in the TL liquid with the dominant contribution of the 2kp oscillating piece. When the long- 
range interaction is present, the correlation decays much more slowly, and thus causes a dominant 
Akp oscillating correlation. This is sometimes referred to as a ID version of the Wigner crystal. 

0) 

We introduce the following interaction in the Hamiltonian (1), 



which smoothly interpolates the short-range interaction (A^O) and the 1/x-type long-range in- 
teraction (A^oo). i) Here the cut-off of x~d is introduced to avoid a singular behavior at short 
distance. The Fourier transformation of (21) gives the modified Bessel function Kq{x), 



Even when the long-range Coulomb interaction is present, we can still diagonalize the Hamiltonian 
as far as the Umklapp and backward scatterings are irrelevant, leading to the formula (9). Note 




(21) 




(22) 



where A =X/d. 
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that Kcip) and Ko{p) now depend on the momentum p, whereas Ks{p) and Kgoip) are independent 
of p, whose value is restricted to Ks{p)=Kso{p)=i as before. 

To discuss the effect of the long-range interaction on the system with the orbital degeneracy, 
we again calculate the single-particle spectrum and the dynamical spin susceptibility. The single- 
particle Green function has the form, 



IJ,=c,o 

p-if 2A+sgn(x+i)^t) if 2A^sgn(x-i;^t) / tt A 



n 
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where 



C^{x,t) = / dp {2A+{p) sm[{x + vi,ip)t}p] + 2A;(p) sin[{x - v^ip)t}p]}, 

Jo P 

(24) 

D^ix,t) = / dp (l + 2np) 

Jo P 

x{2A+{p) (1 - cos[{a; + v^{p)t}p]) + 2A-{p) (1 - cos[{a; - v^{p)t}p])}, 

(25) 

with a'=ao/d and np=l/[exp{v^{p)p/kBT}-i\- In contrast to the short-range case, we cannot 
integrate the charge as well as the orbital parts analytically, so that we evaluate these integrals 
numerically. 

We investigate three cases for which the range of the interaction is long (A'=10^), intermediate 
(A'=1), and short (A'=10~^). As the long-range interaction has no essential effects on the spin 
velocity so it is given by that of the short-range case. We show the momentum-dependence of 
the charge velocity Vc{p) and the TL parameter Kc{p) calculated for the single-orbital case are 
shown in Fig. 5, which are scaled so as to properly reproduce Vc=0.1 and Kc=l/2 in the case of 
the short-range interaction. It is seen that when the long-range interaction is present, they are 
strongly dependent on p in the small p region. 

The single-particle spectrum for the single- as well as two-orbital models is shown in Fig. 6, from 
which we can see how the spectrum changes its profile in the presence of the long-range interaction. 
In the single-orbital case shown in Fig. 6(a), both peaks coming from the charge (ci;~0.01) and 
spin sectors (a;~0.3) are gradually suppressed as the interaction becomes long-ranged (large A ). 
This means that the long-range interaction smoothly drives the system from the TL liquid to the 
ID Wigner crystal. In Fig. 6(b), we show the two-orbital case without the band splitting (A=0). 
Although the spectral function is affected considerably by the orbital degeneracy in the short-range 
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Fig. 5. Charge velocity Vc{p) and TL parameter Kc{p) for the single-orbital model with the long-range interaction. 



interaction, the spectrum loses such characteristic features as the interaction becomes long-ranged. 
Such smearing effect of the spectral function is also seen for the case with the band splitting. The 
single-particle spectrum of the lower and upper bands is shown for A~0.6 in Fig. 7. Therefore, as 
long as the single-particle spectrum is concerned, we would conclude that the long-range interaction 
obscures not only the charge mode but also the other spin and orbital modes in the ID quantum 
liquids. However, this naive conclusion is misleading, which should not correctly capture the effect 
of the long-range interaction. 

To clarify the above point, let us consider the dynamical spin susceptibility. In Fig. 8, we 
show how the long-range interaction affects the dynamical spin susceptibility with and without the 
orbital degeneracy. As is the case for the single-particle spectrum, the peak structure for the charge 
sector is suppressed as the interaction becomes long-ranged. However, it should be noticed that the 
peak structure due to the internal spin and orbital degrees freedom is enhanced as the long-range 
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Fig. 6. Siriglc-particlc spectrum for (a) single-orbital model {vs=tt) and (b) two-orbital model {va=Vo=V3o='K/2). 
The range of the interaction is changed systematically by modifying A'. 



interaction is introduced, as clearly seen in Fig. 8. This characteristic behavior is quite contrasted 
to the results of the single-particle spectrum discussed above. The origin of the enhancement is 
rather clear; i.e. the long-range interaction drives the system to the ID Wigner crystal, which 
almost freezes the charge degrees of freedom and thereby enhances the spin fluctuations effectively. 
When the single-particle spectrum is considered, such enhanced spin fluctuations are again masked 
by the charge excitations which are accompanied by an electron-like excitation. 

Therefore, we come to the conclusion that the internal spin and orbital fluctuations are indeed 
enhanced by the long-range interaction via the formation of the Wigner crystal, although such 
enhancement is obscured in the single-particle spectrum. 

§5. Summary 

We have investigated the effect of the orbital degeneracy and the band splitting on the low-energy 
dynamical properties of ID quantum liquids in a metallic phase. The single-particle spectrum and 
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(a) upper band 




Fig. 7. Single-particle spectrum for the model with short as well as the long-range interaction. The band splitting 
A~0.6, 5^^2.6, 5,0^0.25. 



the dynamical spin susceptibility have been calculated by exploiting the bosonization techniques. 
For the multicomponent TL liquid with the short-range interaction, we have found that the spectral 
function is affected considerably by the orbital degeneracy as well as the band splitting, leading 
to characteristic features of the profile which are quite contrasted to those expected for the Fermi 
liquid. The effects of the long-range interaction have also been discussed. It has been found that 
the charge fluctuations are suppressed while the internal spin and orbital fluctuations are enhanced 
by the long-range interaction via the formation of the Wigner crystal. However, the enhanced spin 
fluctuations are hindered in the single-particle spectrum by the almost frozen charge degrees of 
freedom. It may be interesting to experimentally check whether such enhanced spin fluctuations 
can be indeed observed for ID correlated electron systems, such as the carbon nanotube, etc, for 
which the long-range interaction plays an important role. 
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Fig. 8. Spin susceptibility with the short as well as the long-range interaction. 
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Appendix 

The p-dependence of a{p) and y{p) is obtained as, 

tan[2a(p)] = '^(yip) - y{p)~Y\ 

Vl 

2 _ K,{p)-^ + 1 
^ " Ko' + l ' 
and the TL parameters, Kdp), Vc{p), Kg and Vg, are given by 

K^{p) = [1 + AU{p)ao/TTVi - Vao/TTVi]-^/^ , 
Ko = [l- Vao/TTVi]-'/^ , 
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Vc{p) = [1 + 4C/(p)ao/7rwi - Vao/Trvi]^^^ , 
Vo = wi [1 - Vao/nvif-^^ , 

where we have introduced vi=[vp^ +Vp^]/2 and Av=[vp^-Vp^]/2. Here V is the on-site Coulomb 
interaction. 

Furthermore, the parameters Vi,{p) and K^{p) (i/=c,o) in the diagonahzed Hamiltonian (9) are 
given by, 

X {vc{p)K,{p)Upf + VoKoUpf - 2A«6(P)6(P)} , 
K^{p) = ''c{p)Kcip)Upf + VoKoUp? + 2AvUp)^M 



■«c{v) 



Up? + ±UpY + '^^<i{p)Up) 



^2. . ^ Vc{p)Kc{p)Up? + VoKoiiipf - 2A^6(p)6(y) 
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